Squeezing of optomechanical modes in detuned Fabry-Perot interferometer 
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We carry out analysis of optomechanical system formed by moveable mirror of Fabry-Perot cavity 
pumped by detuned laser. Optical spring arising from detuned pump creates in the system several 
eigen modes which could be treated as high-Q oscillators. Modulation of laser power results in 
parametric modulation of oscillators spring constants thus allowing to squeeze noise in quadratures 
of the modes. Evidence of the squeezing could be found in the light reflected from the cavity. 
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I. INTRODUCTION 

The purpose of gravitational waves detection leads to 
construction of large-scale antennas like LIGO [1, 2], 
VIRGO [3] and GEO [4]. Very high sensitivity of these 
devices is limited by a number of noises. In the low fre- 
quency range ( below ~ 50 Hz) the prevailing sources of 
noise are seismic ones, at middle frequencies (~ 50 — 200 
Hz) thermal noises dominate and in high frequency range 
(over 200 Hz) photon shot noise prevails. However the 
technical improvement of antennas by compensation and 
suppression of these and other noises will allow to achieve 
sensitivity level defined only by quantum noise which for 
continuous position measurement has lowest boundary 
defined by Standard Quantum Limit (SQL) [5-8]. SQL 
is the optimal combination of two noises of quantum na- 
ture: fluctuations of mirror motion caused by random 
photon number falling onto its surface and photon count- 
ing noise. 

One of the ways to overcome the SQL is the imple- 
mentation of so-called optical rigidity (optical spring) 
effect [8-11] which is based on a fact that in a de- 
tuned Fabry-Perot interferometer the circulating power 
and consequently the radiation pressure is dependent on 
the distance between the mirrors. It has been shown in 
a number of papers [12-18] that interferometers using 
optical springs exhibit sensitivity below the SQL. 

In a system utilizing optical spring there are two de- 
grees of freedom (in case of one pump): a mechanical 
one and an optical one. Interaction of these coordinates 
gives birth to several eigen modes (the number of which 
is equal to the system degrees of freedom number) each 
of which is characterized by its own resonance frequency 
and damping. 

Free evolution of the system can be represented as a 
sum of eigen modes each of which is an oscillator with 
its eigen frequency and has corresponding damping. In 
principle one can make transfer from the conventional 
coordinates to eigen ones and consider the evolution of 
the system as evolution of these new oscillators. 

Using the curve of susceptibility of the system one can 
estimate the average energy stored in each equivalent 



'Electronic address: rkh@hbar.phys.msu.ru 



eigen mode oscillator. It is well known [5, 8, 19] that 
an oscillator with average energy exhibits quantum- 
mechanical behavior if the condition 



hujQ 



< 1. 



is met. Here u) and Q are eigen frequency and quality fac- 
tor of oscillator. Estimations [20] show that for system 
such as interferometer Advanced LIGO these conditions 
are well fulfilled. As a consequence one can expect the 
corresponding eigen modes to be observed in quantum 
state despite the fact that the (effective) mass of me- 
chanical oscillator in considered interferometer is equal 
to 10 kg. 

Given this motivation it is interesting to look for any 
experimental scenario giving an ability to bring one of 
the eigen modes to non-classical state. Currently the 
development of techniques aimed on preparation of me- 
chanical resonator in quantum state by means of optome- 
chanical interaction is well underway including investi- 
gations of micromembranes [21], microtoroids [22], op- 
tomechanical crystals [23], pulse-pumped optomechani- 
cal cavities [24] and even large-scale gravitational-wave 
detectors [25]. One of the most important problems in- 
herent to optomechanical devices is relatively large losses 
of mechanical microoscillator. Replacement of material 
spring by optical one may decrease the losses in mechan- 
ical system and thus provide an experimental possibility 
of realisation of quantum regimes. 

For the initial consideration it is also worth to show 
the very possibility of manipulation with optomechan- 
ical modes. To prove this possibility we demonstrate 
the mechanism to perform quadrature squeezing of shot 
noise caused fluctuations in these modes using parametric 
modulation of spring constant which has been previously 
considered as a tool to squeeze fluctuations in cavity mir- 
ror motion [26]. 



II. DESCRIPTION OF MODEL 

We consider a gravitational-wave detector Advanced 
LIGO with a signal recycling mirror (SRM) having 
an amplitude reflectivity p and power recycling mirror 
(PRM); a scheme of antenna is presented in fig. 1. An- 
tenna consists of a Michelson interferometer with addi- 
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FIG. 1: Top: scheme of Advanced LIGO detector. Bot- 
tom: scheme of equivalent Fabry-Perot interferometer with 
one moveable mirror and feedback implemented for stability 
of system. 



tional mirrors forming Fabry-Perot (FP) cavities with 
mean distance L between mirrors in arms. Input mirrors 
have amplitude transmittivities T and masses m and out- 
put mirrors have the same masses m and are completely 
reflective. We assume that all mirrors are lossless. The 
interferometer is pumped by laser having frequency res- 
onant for FP cavities in arms. 

Dynamics of this interferometer could be split into two 
modes: namely differential one and common one. The 
optical fields in the modes represent difference and sum 
of the fields in arms respectively and carry information 
about difference and sum between lengthenings of arm 
cavities. Each mode is characterized by optical detun- 
ing and decay rate introduced by displacement of cor- 
responding recycling mirror: SRM for differential mode 
and PRM for common one. 

It has been demonstrated [15] that each of this modes 
could be equivalently described using a model of FP 
cavity (see fig. 1, bottom). We focus our attention to 
differential mode. The end mirror of equivalent cav- 
ity has mass equal to reduced mass of four mirrors of 
Michelson interferometer /i — ^ and is completely re- 
flective. The displacement of this mirror is equal to 
difference of arms lengthening of initial interferometer: 
y = (x e — Xi) — (y e — y{). Distance separating the mirrors 



is equal to one in initial interferometer, namely L. The 
pump in equivalent scheme should be characterized by 
the same relaxation rate Y and detuning A as in initial 
scheme and mean circulating power P two times bigger 
than one in initial interferometer. 

A detuned pump creates in FP cavity an optical spring, 
that is the radiation-pressure force which depends on the 
distance separating mirrors of the cavity. A system with 
only one optical spring is always unstable because a sin- 
gle pump introduces either negative damping or nega- 
tive rigidity [9-12]. A few ways to avoid instabilities 
have been proposed amongst which there are implemen- 
tation of feedback [15] or utilization of additional pump 
[27, 28]. The latter way has been investigated in details 
and proven experimentally with mirror of gram-scale [29] . 
However in laboratory-scale experiment it should be eas- 
ier to utilize a proper feedback. Hence we consider a 
scheme of FP cavity with one pump and a feedback in- 
stead of two pumps. 

Input optical fluctuations are described by annihila- 
tion operator a 1 , operator a° describes output fluctua- 
tions. To avoid instability caused by optical spring we 
measure phase quadrature aJJ of output fluctuations with 
homodyne detector and send its derivative over time as 
a feedback force to the moveable mirror. 

The dynamics of interferometer could be described in 
terms of two degrees of freedom. One of those is annihi- 
lation operator b of optical fluctuations inside the inter- 
ferometer (actually, its amplitude quadrature b\) and an- 
other is mechanical displacement y. Equations of motion 
for this system could be written as follows (in capacity 
of mechanical coordinate we use dimensionless displace- 
ment z): 



h+gh + 26a + Az = v x ; 
z — Ab\ + ab\ = V2- 



(la) 
(lb) 



Here we use dimensionless parameters defined the same 
way as in [20, 28] . Parameter g stands for optical relax- 
ation rate, coupling coefficient A is proportional to pump 
power, a is the coefficient of feedback. The right parts of 
equations represent fluctuational forces acting on corre- 
sponding degree of freedom. As of optical one the force 
V\ represents input vacuum fluctuations and for mechan- 
ical one vi describes vacuum fluctuations reflected by the 
cavity and transmitted by feedback. We neglect all other 
fluctuational forces including thermal and seismic ones 
(although this is only an assumption for initial consider- 
ation it is quite reasonable for Advanced LIGO [2]). For 
exact definitions see Appendix A. 

The equations (1) look similar to ones for oscillator 
coupled to a free mass (or two coupled oscillators one of 
which has partial frequency equal to zero). However we 
would note the difference in signs in front of the coupling 
terms (proportional to A): in the case of conventional 
coupled oscillators these signs coincide; in our case the 
signs differ. This is the reason of instability which ap- 
pears in absence of feedback and this is also the reason 
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of possibility for eigen frequencies of such system to coin- 
cide (so called double resonance case [13, 16, 17] — let us 
remind that in case of two coupled oscillators coincidence 
of eigen frequencies is impossible). 



III. EIGEN MODES DECOMPOSITION 

As in the case of conventional coupled oscillators one 
can treat the system in terms of eigen modes. The evo- 
lution of oscillator system could be written as a sum of 
eigen oscillations: 



Hi is unnecessary). Vectors W{ are defined as follows: 



-lUJit 



(2) 



Index in the sum runs from 1 to the number of system 
degrees of freedom (2 for set (1)). This equation could be 
treated as transformation from conventional coordinates 
(b\; z) to the eigen ones. As in eqn. (2) one can present 
evolution of each of eigen modes as combination of os- 
cillations at corresponding eigen frequencies uji and slow 
(compared to these oscillations) changing of amplitude. 
In case of free evolution (zero right parts in eqns. (1)) 



quantities <?j 



-fit ; 



i.e. they freely decay with relax- 



ation rates ji into the equilibrium values. Vectors v*i rep- 
resent distribution of amplitudes in corresponding eigen 
modes. 

Suppose that we have a possibility to modulate the 
power of pumping laser. This will result in modulation 
of coupling coefficient and hence it will cause shifting 
of eigen frequencies of system. It is well known [30] 
that modulation of eigen frequency of oscillator with fre- 
quency twice bigger than its own one results in squeezing 
of noise in its quadratures. In this paper we show that 
modulation of the pumping power performs squeezing of 
noise in quadratures of eigen modes amplitudes. 

Consider we harmonically modulate the pumping 
power at frequency 2p so coupling coefficient in (1) de- 
pends on time: A — > A(l + 2\m\ cos(2pt + </>)). Assum- 
ing the quantities gi to be slow (i.e. not to significantly 
change on times compared to mechanical periods 2ir/uii) 
one can plug the expression (2) into system (1) and to 
get rid of rapidly oscillating terms by averaging over a 
period 2tt/p. This procedure results in shortened equa- 
tions for amplitudes gi. For one with number j (J = 1, 2) 
the equation takes form: 



ujj—2p)t—i(p 



feOe*"*. (3) 



Here we have introduced a set of vectors Ilj built to be or- 
thogonal to all Vi except one: (Tljili) = Sji, where 5ji is a 
Kronecker delta (in the case of two conventional coupled 
oscillators vectors Vi are itself orthogonal and system of 



A 



Underlining used for is to emphasise that this quan- 
tity is obtained in right part of equation by procedure 
of keeping only slow (in comparison to terms oscillating 
with frequency p) quantities. To make Vjt luljt fulfilling 
this criterion one should keep in Vj only spectral compo- 
nents close to uij . 

From now on we will focus on the amplitude of one of 
modes taking into account that all the results obtained 
for this mode are similar for other ones. For clarity let 
us set j = 1 thus considering the first mode. Also we 
suppose the modulation to happen at frequency twice 
bigger than the one of first mode: p = u)\ . 

First let us consider the simplest case when difference 
between eigen frequencies u>\ — u>2 is large in respect to 
decay rates 7$ . In this case the exponential multipliers in 
eqn. (3) should be considered as fast oscillating ones and 
the equations for amplitudes gi decouple. In this case 
the equation for g\ takes the following form in spectral 
domain 

[71 ~ ix]gx(x) + engK-x) = = h(x). (4) 

Here x stands for normalized spectral frequency (see ap- 
pendix A), 71 is the decay rate of first mode; en is a quan- 
tity proportional to modulation strength, the general ex- 
pression for it could be easily deduced from eqn. (3): 

-i\m\(U* j w i )e i ' / ' 
e * = " -to, • 

We also write the equation for g\(—x) by hermitian 
conjugation of eqn. 4 and replacement x — > —x. Taking 
proper combinations of these equations yields the expres- 
sions for corresponding quadratures of g\ . For simplicity 
let en = e n = |en| which is achievable by proper choos- 
ing </>, in this case the mentioned combinations reduce to 
sum or difference: 



{9i+9[) 

(.91 -g\) 



7i 
7i 



I en I ~ ix 
kill - ix 



= h + fh 
= h + fb 



It is obvious from these equations that sum quadrature 
is squeezed due to parametric modulation and the dif- 
ference quadrature is antisqueezed. It is also could be 
seen that quantity |en| has sense of additional damp- 
ing (positive or negative for different quadratures) in- 
troduced by modulation. In absence of modulation en 
should be equal zero and there should be no difference 
between quadratures. 

One can also estimate the critical level of modulation 
characterized by value m c , this is a level when negative 
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damping added to the differential quadrature gi — g\ be- 
comes equal to its own damping: 71 — |en| = thus mak- 
ing quadrature instable. Given this value of modulation 
coefficient damping in the sum quadrature is twice bigger 
than its own damping 71. This means that quadrature 
noise squeezing is limited by factor of two. 

In general case of arbitrary eigen frequencies cli 12 the 
equations (3) could not be solved in such an obvious way, 
but in this case one can rewrite these equations to switch 
from spectral components gj (x) to the same components 
but shifted by frequency, namely gj(x — (ojj — wi)). The 
equations take form 

gj(x — (uij - u>x)) ■ [ — ix + i(u>j - u>x) + jj] + 

+ ^9i(-x-(ui-"i)V j i = -^(ftjEjip + x)). (5) 



This set represents a system of linear algebraic equa- 
tions and could be solved by applying Kramer's rule. Us- 
ing the solution we obtain the expression for plus or mi- 
nus quadratures of j-th mode Gj defined as follows 
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FIG. 2: Spectral density of output field quadratures A\ 
(defined in eqn. (6)). Trace 1: spectral density of A[ + \x) 
and A[ \x) in absence of parametric excitation (m = 0), 
plots coincide. Traces 2 and 3: spectral density of squeezed 
quadrature A[ + \x) for modulation coefficient equal to m = 
0.35m c and m — 0.7m c . Traces 4 and 5: spectral density 
of anti-squeezed quadrature A[~\x) for the same values of 
modulation coefficient. 



Gf\x) 



9j(x)±9j(-x) 
V2 



and calculate spectral densities of these quadratures. 

In coincidence with expectations in absence of para- 
metric modulation both plus and minus quadratures of 
first mode have equal spectral densities. If modulation is 
applied the plus quadrature G^ appears to be squeezed 
and the minus quadrature G\ ' to be antisqueezed. The 
bigger is modulation coefficient the more significant the 
(anti-) squeezing effect is. It could be easily shown that 
squeezing is also limited by a factor of two in this case. 



frequency ujj with proper phase: 
a(t) cos(ujjt + 4>m) — > a(x + Wj)e l * M + a(x — Wj)e~ 



i<t>t. 



Setting 



= or 



>M 
l(±) 



= ? we obtain Ay' or A) 



Estimations of Aj spectral density demonstrate that 
it is possible to see squeezing in output light. The spec- 
tral densities of these quantities are plotted in fig. 2 be- 
ing normalized to the shot noise level. The limitation of 
squeezing by factor of two shows up in output light too. 



IV. READOUT 

The important question is whether it is possible to see 
this squeezing in output light. The answer to this ques- 
tion is positive. 

Output fluctuations a° are defined by input fluctua- 
tions a 1 and light inside cavity b, in particular for phase 
quadrature of a° the following expression is valid: 

a° 2 = -a\ - 2VIV&1. 

Plugging the solution obtained in previous section for 
gx,2 into the equation (2) we can thus express the output 
fluctuations through quantities gi. It is easy to show (see 
appendix B) that quantity 

is proportional to the quadrature G^ ± \x). In time do- 
main measurement of this quantity is equivalent to mul- 
tiplication of measured phase quadrature with cosine of 



V. CONCLUSION 

In this paper we have shown the possibility of quadra- 
ture noise squeezing of eigen modes amplitudes in an op- 
tomechanical system. As a model of latter we have used a 
Fabry-Perot cavity with moveable mirror and a detuned 
pump. The eigen modes in this system are defined by in- 
teraction between optical and mechanical degrees of free- 
dom. 

Detuned pump transforms the mirror which is initially 
a free-mass into a harmonic oscillator with spring coeffi- 
cient provided by means of optical rigidity. In comparison 
with material rigidity peculiar to usual optomechanical 
devices such as microtoroids or membranes etc [21-24] 
optical rigidity is characterized by very low noises level. 
The another advantage of using optical springs is the 
very ability of spring constant manipulation by means 
of pump power modulation. The fact that main noise 
source in planned Advanced LIGO is the quantum one 
gives us an ability to talk about quantum noise squeezing 
in quadratures. 
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Despite the fact that we consider a feedback in our 
model, the feedback itself is unnecessary to achieve the 
squeezing as it serves only to stabilize the eigen frequen- 
cies of the system. The stabilization could be performed 
by means of an auxiliary pump with power significantly 
lower than one in main pump. However implementation 
of feedback is more feasible in small-scale experimental 
setups like gravitational-wave antenna prototypes [31- 
33]. 

Realization of proposed scheme in experiment could 
demonstrate non-classical behavior of optomechanical 
mode which could be useful for tests of quantum mechan- 
ics applied to macroscopic mechanical objects. We would 
like to underline it is not pure optical or pure mechani- 
cal, but an optomechanical degree of freedom that may 
exhibit quantum behavior. Also a quantum state if cre- 
ated in an eigen mode could be transmitted into another 
one (for example, by modulation of pump with difference 
frequency {u>2 — ^>i))- This gives us a potential play- 
ground for creation and transmission of quantum states 
between optomechanical modes which could be useful for 
problems of quantum information. 
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as follows. 



V2n 



2V2T 



~' w 2hr m x ' ,/r 2 + A 2 ' y Vr2 + A2' 



.4 = 



Vr 2 + A 2 V v Lc 



2^2 /2/zAT 
ct = — ath \ / — - — ■ 

Vr 2 + A 2 V h 

Switching to dimensionless parameters from dimensional 
ones using the definitions listed above and from frequency 
domain to time one following rule — iQ — > d t yields sys- 
tem (1). 

For numerical calculations through this paper we have 
used the following set of parameters 

A = 0.90; g = 0.009; a = -0.01. 



Appendix B: Readout 

In output field we measure phase quadrature aSJ which 
is proportional to b\ that carries information about 
squeezing. Further we discuss which combination of spec- 
tral components of b\ contains squeezing evidence taking 
into account that corresponding combination of aJJ com- 
ponents represents the same combination with some in- 
put fluctuations a\ added. 

For spectral components of b\ shifted in frequency do- 
main by amount of uij one can write the following expres- 
sion using eqn. (2): 

bi(x+Uj) = ^TVi(gi(x- [ui-<jj])+gl(-x- [wi+Uj])). 



Appendix A: Details of notations in initial equations 



Below we write down the system of equations of motion 
for interaction of light quadrature and mirror displace- 
ment in spectral domain using dimensional parameters. 



Pujo 
hi 2 



A 



[ai(r-iO) + a|A] ; 



-O 2 j/ \J Phuiobi — — i51o!fb(— a\ — 2%/Tt&i). 



Here f2 is spectral frequency, r is the time it takes light 
to travel between mirrors (one way), afb is a coefficient 
of feedback, i.e. the force fed to the mirror is equal to 
/fb = afb^^- The dimensionless parameters are defined 



Here we use freedom in definition of Vi and let Vi — 
(^i)i = (^i)i- An equation similar to one above could 
be written for b\ (x — ujj): 

h (x-u)j)=^2 V i (dl [-x-[u>i-Uj])+ 9i (x -[ui + uij])). 



We now take sum or difference of these quantities and 
obtain expression containing sum of quadratures Gj 
that are of our interest 



B 



(±) _ h(x + ujj) ±bi(x — uij) 



V2 



+ (gi(x - (wi + Uj)) ± g\{-x - (w s + Uj))) = 

i 

Remind that we assumed gi to be slow amplitudes 
which means that spectral components of these quan- 
tities are situated close to zero frequency. Note that due 
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to resonant multiplyers (expressions in square brackets) 
in eqn. (5) one can estimate the width of band contain- 
ing §i spectra with order of corresponding mode decay 
rate 7$ which is much smaller than difference between 
eigen frequencies u>i — u)j (for i ^ j). This considera- 
tion is also valid for quadratures Gi hence in definition 
of all summands differ from zero in different fre- 

quency bands (not overlapping because of narrowness), 
or in other words in each frequency band there is not 
more than one non-zero summand in this definition. In 
particular at frequencies close to zero one can use quite 
exact expression 



V jG f 



Or). 



As we actually measure a 2 we need to consider a com- 
bination of spectral components of this quantity contain- 
ing Bj. It is obvious that corresponding combination 
is given by the following expression 



(±) _ a° 2 {x + uj )±a 2 {x-uj j ) 



a l 2 {x + LOj) ± a 2 (x — LOj) 

71 



2^TT-B<f\x) 



which represents a sum of desired quadrature and input 
fluctuations. 
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